0–π transition induced by the barrier strength in spin superconductor Josephson junctions*

Project supported by the National Key R&D Program of China (Grant No. 2017YFA0303203) and the National Natural Science Foundation of China (Grant No. 11474149).

Zeng Wen1, Shen Rui1, 2, †
National Laboratory of Solid State Microstructures and School of Physics, Nanjing University, Nanjing 210093, China
Collaborative Innovation Center of Advanced Microstructures, Nanjing University, Nanjing 210093, China

 

† Corresponding author. E-mail: shen@nju.edu.cn

Project supported by the National Key R&D Program of China (Grant No. 2017YFA0303203) and the National Natural Science Foundation of China (Grant No. 11474149).

Abstract

The Andreev-like levels and the free energy of the spin superconductor/insulator/spin superconductor junction are obtained by using the Bogoliubov–de Gennes equation. The phase dependence of the spin supercurrents exhibits a 0–π transition by changing the barrier strength. The dependences of the critical current on the barrier strength and the temperature are also presented.

1. Introduction

In the conventional superconductor, the electrons are condensed into the Cooper pairs, forming a Bardeen–Cooper–Schrieffer (BCS) ground state.[1] Each Cooper pair carries two units of the electron charge and no spin polarization for spin-singlet pairing. Being the counterpart of the charge superconductor, there also exists a quantum state, called a spin superconductor (SSC), where the condensed pair is composed of an electron and a hole rather than two electrons.[25] The electronlike spin-up carriers and the holelike spin-up carriers attract each other due to the Coulomb interaction[6] and form the Cooper pairs in the SSC, which are the charge neutral spin-triplet pairs carrying two units of the electron spin. About four decades ago,[7,8] electron–hole Cooper pairs were predicted in a double-layer system, where the electron in one layer and the hole in the other layer are separated by an insulator. The insulating layer is proposed to be thin enough to allow for the strong Coulomb interaction between the electrons and the holes leading to an effective attraction but thick enough to prevent the electron–hole recombination. Experimentally, the existence of the exciton condensate in the semiconductor double quantum well systems was also reported.[9,10] In addition, the spin superconductivity can exist in various systems, such as the spin-polarized triplet exciton systems of graphene,[2,3] Bose–Einstein condensate of magnetic atoms,[1113] Bose–Einstein condensate of magnons and spinons,[1416] and so on. Very recently, the signature of spin superconductivity has also been observed experimentally in canted antiferromagnet Cr2O3.[17]

The conventional Josephson effect is the phenomenon of the supercurrent across a device known as the Josephson junction, which consists of two superconductors coupled by a weak link.[1823] The charge supercurrent across the junction is mediated by the Andreev bound states formed by the electrons and holes inside the normal segment of the junction.[2431] As an analogy of the conventional Josephson junction, the SSC Josephson junction can be achieved by a weak link between two SSCs. There are both the spin-conserved and the spin-flip reflections at the normal-metal/SSC interface,[3] which are the counterparts of the normal reflection and the Andreev reflection at the normal-metal/superconductor interface, respectively. The repeated spin-flip reflections in the SSC Josephson junction, just like the Andreev reflection in the conventional Josephson junction, will lead to the formation of Andreev-like bound states and a dissipationless pure spin current.[32]

In this paper, we investigate an SSC Josephson junction by employing the Bogoliubov–de Gennes (BdG) approach, which is completely in parallel to the treatment in the conventional Josephson junction. The Andreev-like levels and the free energy of the junction are obtained. It is shown that there is a pure spin supercurrent driven by the phase difference of the junction. The current phase relation can exhibit a 0–π transition by tuning the barrier strength. We note that there are two differences between the SSC Josephson junction and the conventional one. First, the supercurrent in the SSC Josephson junction is a pure spin current rather than a charge current. Second, the 0–π transition is tuned by the barrier strength rather than the magnetic middle layer.[3336] The dependences of the critical current on the barrier strength and the temperature are also discussed.

The rest of this paper is organized as follows. The formalism is presented in Section 2. The numerical results are discussed in Section 3. Finally, the main results are summarized in Section 4.

2. Formalism

The model under consideration is a one-dimensional junction oriented along the x axis. The junction consists of two SSCs sandwiched by an insulating layer at x = 0. The SSC can be described by a minimal two-band model,[6] where an electron in the spin-up conduction band and a hole in the spin-down valence band can be condensed into a Cooper pair due to the Coulomb interaction. Both the electrons in the spin-up band and the holes in the spin-down band possess the spin of +ħ/2, thus the Cooper pair is charge neutral and has a spin of ħ. The dispersion of the spin-up conduction band and that of the spin-down valence band are assumed as ε+ = p2/2mM0 and ε = −p2/2m+M0, respectively, with p being the momentum, m the mass of the electron, and M0 the exchange energy. The spin-down conduction band and the spin-up valence band are far from the Fermi level and are irrelevant to the transport.

In the basis of (ψ+, ψ)T, with ψ± being the field operators of the electrons in band ε±, the BdG Hamiltonian of the SSC reads

where Δ is the pair potential between the electrons in band ε+ and the holes in band ε. The pair potentials in the left and right SSCs have a phase difference φ. The insulating layer is modeled by a barrier potential (x). The eigenfunction of the BdG Hamiltonian can be obtained as
where
with E being the energy of the quasiparticles. The wave vector is written as . For the states with E > |Δ|, the wave vectors are real, while for the subgap states, the wave vectors have a small imaginary part so that the wave functions properly decay at the infinity. The coefficients a, b, c, and d can be determined by matching the boundary conditions
where σz is the Pauli matrix.

The Pauli matrix in Eq. (3) indicates that the boundary condition in an SSC junction is quite different from the one in a conventional superconductor junction. For the conventional superconductor, the two diagonal blocks of the BdG Hamiltonian describe the electrons and the holes, respectively. Both the kinetic energy and the barrier potential possessed by the electrons are opposite to those by the holes so that there is no Pauli matrix in the boundary condition for a nonmagnetic junction. For the SSC, the two diagonal blocks of the BdG Hamiltonian describe the electrons in the conduction band and the electrons in the valence band, respectively. The conduction electrons have opposite kinetic energy to the valence electrons but they feel the same barrier potentials, which results in the Pauli matrix in the boundary condition (3). Therefore, the barrier potential in the SSC junction may play a role like the magnetic scattering in the conventional superconductor junction, and may induce a 0–π transition in the SSC Josephson junction as discussed below.

Boundary conditions (3) are the linear equations of the coefficients a, b, c, and d. The nontrivial solution demands a secular equation, which in turn gives four Andreev-like levels in the junction as

where Δ0 is the magnitude of the pair potential and z = 2mV/ħ2k0 is the dimensionless barrier strength with . The φ-dependent part of the free energy of the junction is given by[33]
where β = 1/kBT with temperature T and Boltzmann constant kB. The dissipationless current carried by the Cooper pair in the junction is determined by[26]
where f(εi) is the Fermi–Dirac distribution function. Since the Cooper pair in the SSC is charge neutral and has a spin of ħ, the current in Eq. (6) is a pure spin supercurrent.

3. Results and discussion

The free energy of the junction is plotted in Fig. 1, which is a function of the phase difference and the barrier strength. When the minimum of the free energy is at φ = 0 or φ = π, the junction is called a 0-junction or a π-junction, respectively. One finds from Fig. 1(a) that the junction is a 0-junction with the small barrier and a π-junction with the large barrier, resulting in a 0–π transition induced by the barrier strength. This feature is different from the conventional Josephson junction, where the 0–π transition is induced by the magnetic scattering.[37]

Fig. 1. (color online) (a) Free energy as a function of φ at T = 0.2Tc, where Tc is the superconducting critical temperature. The red solid line and the purple dash-dotted line are for z = 0.3 and z = 4, respectively. The blue dashed line and the orange dotted line are for z = 1.5 and z = 1.7, respectively. (b) Free energy as a function of z. The red solid line and the blue dashed line are for φ = 0 and φ = π, respectively.

The 0–π transition point zt can be found as the cross point of the curves F(z,φ = 0) and F(z,φ = π) as shown in Fig. 1(b). For the barrier strength far from the transmission point, the free energy has only two extrema. The minimum point is at φ = 0 and the maximum point is at φ = π or vice versa, as shown by the red solid line and the purple dash-dotted line in Fig. 1(a). With increasing barrier strength, the free energy at φ = 0 increases and the one at φ = π decreases, until they are equal at z = zt. During this process, there must be a new extreme point present in the region 0 < φ < π to prevent the free energy from being a flatline at z = zt, as shown by the blue dashed line and the orange dotted line in Fig. 1(a), and the two curves of the free energy with the barrier strength near the transition point.

The 0–π transition can also be indicated by the phase dependence of the Josephson spin current, which is depicted in Figs. 2(a) and 2(b). For the 0-junction with z = 0.3, the phase dependence of the Josephson spin current is a sinusoidal function with the maximum spin current occurring in the region of 0 < φ < π. For the π-junction with z = 4, the current–phase relation is a negative sinusoidal function with the maximum spin current occurring in the region of π < φ < 2π. For the junctions with the barrier strength near the transition point zt, the current–phase relation deviates from the sinusoidal form and has additional zero points besides φ = 0 and φ = π, which is a result of the extra extreme point in the free energy. However, the current–phase relation is still a 2π periodic function of φ. The maximum spin current still occurs in the region of 0 < φ < π for the 0-junction and in the region of π < φ < 2φ for the π-junction, as shown by the red solid curve and the blue dashed one in Fig. 2(b), respectively.

Fig. 2. (color online) The phase dependence of the Josephson spin current. (a) The red solid line and the blue dashed line are for z = 0.3 and z = 4, respectively. (b) The red solid line and the blue dashed line are for z = 1.5 and z = 1.7, respectively.

The critical spin current Ic, which is defined as the maximum Josephson spin current, is shown in Fig. 3 as a function of the barrier strength. One finds that there is a cuspidal point in the z dependence of Ic, which also indicates the 0–π transition. The current acquires its maximum value at the critical phase difference φ0, which is a function of z. Therefore, we have Ic = I(φ0(z),z) with its derivative written as

Since φ0 is the extreme point of the current, we have (∂I/∂φ)|φ0 = 0 and dIc/dz = (∂I/∂z)|φ0. For the 0-junction and the π-junction, φ0 is less than π and greater than π, respectively, as shown in Fig. 2. Therefore, there is always a jump of φ0 and a discontinuity in the derivative of Ic around the transition point zt, as shown in Fig. 3. Such a discontinuity is also shown in the 0–π transition in the conventional superconductor/ferromagnet/superconductor junctions, however, the transition is driven by the thickness of the junction there.[38]

Fig. 3. (color online) The critical current as a function of z.

Finally, we present the temperature dependence of Ic in Fig. 4 by assuming the BCS-like temperature dependence of the pair potential. One finds that the critical spin current decreases naturally with temperature increasing and vanishes at the superconducting critical temperature.

Fig. 4. (color online) The temperature dependence of Ic. The red solid line, the blue dashed line, the orange dotted line, and the purple dash-dotted line are for z = 0.3, z = 1.5, z = 1.7, and z = 4, respectively.
4. Summary

We investigate the dc Josephson effect in an SSC Josephson junction. In contrast to the conventional Josephson junction, the Josephson current is a pure spin current and the 0–π transition is driven by the barrier strength rather than the magnetic mechanism. The signals of the 0–π transition are shown in the free energy, in the current-phase relation, and in the critical Josephson current.

Reference
[1] Bardeen J Cooper L N Schrieffer J R 1957 Phys. Rev. 108 1175
[2] Sun Q F Jiang Z T Yue Y Xie X C 2011 Phys. Rev. 84 214501
[3] Liu H W Jiang H Xie X C Sun Q F 2012 Phys. Rev. 86 085441
[4] Sun Q F Xie X C 2013 Phys. Rev. 87 245427
[5] Bao Z Q Xie X C Sun Q F 2013 Nature Communications 4 2951
[6] Lv P Guo A M Li H Y Liu C X Xie X C Sun Q F 2017 Phys. Rev. 95 104516
[7] Shevchenko S I 1976 Fiz. Nizk. Temp. 2 251
[8] Lozovik Y E Yudson V I 1975 JETP Lett. 22 274
[9] Sivan U Solomon P M Shtrikman H 1992 Phys. Rev. Lett. 68 1196
[10] Butov L V Zrenner A Abstreiter G Bohm G Weimann G 1994 Phys. Rev. Lett. 73 304
[11] Chang M S Qin Q S Zhang W X You L Chapman M S 2005 Nat. Phys. 1 111
[12] Lin Y J Garci K J Spielman I B 2011 Nature 471 83
[13] Hamner C Zhang Y P Khamehchi M A Davis M J Engels P 2015 Phys. Rev. Lett. 114 070401
[14] Giamarchi T Rsegg R Tchernyshyov O 2008 Nat. Phys. 4 198
[15] Mazurenko V V Valentyuk M V Stern R Tsirlin A A 2014 Phys. Rev. Lett. 112 107202
[16] Kimura S Kakihata K Sawada Y Watanabe K Matsumoto M Hagiwara M Tanaka H 2016 Nat. Commun. 7 12822
[17] Yuan W Zhu Q Su T Yao Y Y Xing W Y Chen Y Y Ma Y Lin X Shi J Shindou R Xie X C Han W 2018 Science Advances 4 eaat1098
[18] Josephson B D 1974 Rev. Mod. Phys. 46 251
[19] Josephson B D 1962 Phys. Lett. 1 25
[20] Anderson P W Rowell J M 1963 Phys. Rev. Lett. 10 230
[21] Zhai J Q Li Y C Shi J X Zhou Y Li X H Xu W W Sun G Z Wu P H 2015 Chin. Phys. Lett. 32 47402
[22] Wu B H Feng X Y Wang C Xu X F Wang C R 2016 Chin. Phys. Lett. 33 17401
[23] Su F F Liu W Y Xu H K Deng H Li Z Y Tian Y Zhu X B Zheng D N Lv L Zhao S P 2017 Chin. Phys. 26 060308
[24] Beenakker C W J 1991 Phys. Rev. Lett. 67 3836
[25] Beenakker C W J Houten H V 1991 Phys. Rev. Lett. 66 3056
[26] Beenakker C W J Houten H V 1991 Proceedings of the International Symposium on Nanostructures and Mesoscopic Systems Santa Fe Kirk W P 20 24
[27] Giuliano D Affleck I 2013 Journal of Statistical Mechanics 02 02034
[28] Wang Z Y Shen R 2010 Chin. Phys. 19 087401
[29] Hamidreza E 2014 Chin. Phys. 23 057402
[30] Li X W Liu D Bao Y H 2013 Chin. Phys. 22 127401
[31] Wang P Xie W Hu L Liu X Zhao X J He M Ji L Zhang X Yan S L 2013 Chin. Phys. 22 057402
[32] Wang B G Peng J Xing D Y Wang J 2005 Phys. Rev. Lett. 95 086608
[33] Annunziata G Enoksen H Linder J Cuoco M Noce C Sudb A 2011 Phys. Rev. 83 144520
[34] Lin Y J Garci K J Spielman I B 2011 Nature 471 83
[35] Meng Q L Shivamoggi V Hughes T L Gilbert M J Vishveshwara S 2012 Phys. Rev. 86 165110
[36] Kwon H J Sengupta K Yakovenko V M 2004 Euro. Phys. J. 37 349
[37] Ryazanov V V Oboznov V A Rusanov A Y Veretennikov A V Golubov A A Aarts J 2001 Phys. Rev. Lett. 86 2427
[38] Oboznov V A Bolginov V V Feofanov A K Ryazanov V V Buzdin A I 2006 Phys. Rev. Lett. 96 197003